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CALIBRATION AND ERROR CORRECTION IN MULTI-CHANNEL IMAGING 

BACKGROUND OF THE INVENTION 

5 Field of the Invention 

[0001] The present invention relates to means for calibrating 
multi-channel imaging systems. In particular, the invention 
provides an approach for measuring and correcting the field- 
dependent image distortion, attenuation, and polarization rotation 
10 that may be present in multi-channel imaging systems to improve the 
accuracy of interf erometric measurements of optical wavefronts. 

Description of the Related Art 

[0002] Multi-channel imaging systems can be used for the 
15 measurement of spectroscopic, polarimetric, and/or interf erometric 
properties of objects by simultaneously acquiring a plurality of 
images on either single or multiple detector arrays. See, for 
example, U.S. Patents No. 5,926,283, No. 5,982,497, No. 4,575,248, 
- - No. 5, 589, 938, No. 5, 663, 793, . No ._ 5, 777, 741, No. 5, 883, 717, No. 
20 4, 624, 569, and No. 6,304, 330. Data from the multiple images can be 
used for qualitative comparison or can be combined in quantitative 
ways. In the case of polarimetry and interf erometry, the data at 
each corresponding image pixel are combined through a mathematical. 
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algorithm to solve for a desired fundamental optical parameter 
(e.g., the optical path difference or the polarization state) . The 
registration and losses between images can be extremely important 
to the accuracy of the calculation. 

5 [0003] For ideal imaging systems, the registration between images 
can be accomplished through simple lateral shifting of the images 
digitized from the detector arrays. Thus, the registration can be 
expressed as a simple x and y offset for each image. However, even 
highly corrected real-world optical systems will contain field- 

10 dependent image distortion, which may not be perceptible to the 
human eye but can cause significant measurement errors. To 
mitigate these errors, it is first necessary to measure the image 
distortion and then to construct an algorithm to adjust the images 
and correct for the distortion. The prior art teaches methods for 

15 transforming images into alternate coordinate systems for computer 
graphics displays, but it is silent with respect to applying these 
algorithms to multi-channel imaging systems to improve measurement 
accuracy. See U.S. Patent No. 6,249,289; L. G. Brown, "A Survey of 
Image Registration Techniques," Columbia University, Computer 

20 Science Department Publication, New York, NY, 1991; and R. C. 
Gonzalez et al., Digital Image Processing , Addison-Wesley, Reading, 
MA, 1987. 
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[0004] Another problem that can greatly affect the accuracy of 
multi-channel imaging systems is the non-linear detector response 
between corresponding image pixels. This can be caused by 
electrically dead pixels on the detector array, non-linear electro- 
optical response, or obscuration due to contaminants such as dust 
or coating defects within the optical imaging system. The prior 
art describes methods for identifying phase-calculation algorithms 
that are robust with respect to non-linear detector response; 
however, these algorithms assume the same detector pixel is used 
for each phase-shifted image (see J. Schmit et al., "Extended 
averaging technique for derivation of error-compensating algorithms 
in phase-shifting interf erometry, " Applied Optics, Vol. 34, p. 
3610, July 1995) . The. prior art describes methods for calculating 
data for bad or missing pixels by using nearest neighboring pixels; 
however, these also assume the same detector pixel is used for each 
phase-shifted image (see C. K. Hong et al., "Least-squares fitting 
of the phase map obtained in phase-shifting .electronic speckle 
pattern interf erometry, " Optics Letters, Vol. 20, p. 931, April 
1995) . Further, these algorithms reduce spatial resolution. In 
addition, prior-art "methods for calibrating multi-channel 
interferometer systems describe methods for aligning the optical 
system to minimize registration errors between images but do not. 
account for geometrical image distortion or pixel dependent non- 
linearities that may be present in the system, which can 



significantly limit accuracy. See T. D. Upton et al., "Optical and 
electronic design of a calibrated multichannel electronic 
interferometer for quantitative flow visualization," Appl . Opt.. 
Vol. 34, No. 25, 1995; and C. Koliopoulos et al., "Simultaneous 
phase shift interferometer," SPIE Vol. 1531, pp. 119-127, 1991. 

[0005] Schwider et al. (Applied Optics, Vol. 22, pp. 3421-3432, 
1983) and Schmit et. al. (Applied Optics Vol. 34, pp. 3610-3619, 
1995) describe the use of extended averaging by introducing a know 
phase-shift between measurements to eliminate residual phase- 
dependent systematic error in conventional temporal phase-shifting 
interferometers." Multi-channel interferometers are typically 
employed in situations where the relative phase between the test 
and reference arms is unstable (e.g., due to vibration or 
turbulence) . Therefore, introduction of a precise phase-shift 
between measurements is not possible, in general, and these methods 
are not adequate. 

[0006] Schwider et. al. (Applied Optics, Vol 28, No. 18, pp. 3889- 
3892, 1989; and Applied Optics, Vol. 38, No. 4, pp. 655-659, 1999) 
also describe an a posteriori technique for calibrating errors in 
phase-shifting interferometers. This method requires the 

introduction of several tilt fringes, fitting the measured data to 
determine the phase-dependent error and then reprocessing the data 
to . remove the error. In general, this method requires 



recalibration for every new tilt orientation and does not work for, 
optical alignments near null-fringe condition where interferometers 
are typically operated to minimize wavefront error. Thus, a 
general method for the elimination of residual systematic errors in 
multi-channel interferometer systems is still needed and very 
desirable. 



BRIEF SUMMARY OF THE INVENTION 



[0007] In view of the aforementioned limitations of the prior art, 
it is one object of the present invention to provide a means for 
5 measuring the field-dependent image distortion in a multi-channel 
imaging system and a means for correcting subsequently measured 
images. 

[0008] A further object of the invention is to provide a means for 
storing the image distortion information in a compact file format 
10 that can be rapidly recalled and applied to subsequent 
measurements . 

[0009] A further object of the invention is to provide a means for 
detecting dead or obscured pixels in one image and approximating 
the value of such a pixel through the use of. values measured at 
15 corresponding pixels in the other images. 

[0010] A further object of the invention is to provide a means for 
measuring the relative photo-response of corresponding pixels in 
each image/" storing" coefficients, and applying a correction for 
normalizing the response of each pixel. 

20 [0011] An additional objective of the invention* is to provide a 
means for reducing phase-dependent systematic measurement error of 
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interf erometric phase measurements by averaging multiple 
measurements with random phase offsets of the reference path 
length. 



[0012] In accordance with these objectives, the invention 
5 consists of measuring field-dependent geometric distortion for 
each image produced by the multi-channel system and calculating 
correction coefficients to remove the distortion from the 
measurement. According to one aspect of the invention, geometric 
distortion , is determined by projecting a known pattern into the . 

10 multi-channel imaging system and measuring the position of the. 

pattern in each image with respect to a predetermined reference. 
Such measurement is accomplished by parsing each image into a 
. plurality of regions with predetermined spatial coordinates, 
finding the absolute or relative offset of the known pattern in 

15 each region, and "performing a polynomial surface fit of the 

pattern iri each region to enable interpolation between region 
patterns. The polynomial coefficients are then stored and 
subsequently used to correct the measurement data. According to 
another approach, geometric distortion is measured using cross- 

20 correlation of the images with an arbitrary or a random 

calibration test pattern. This method only enables the relative 
correction between images, but it has the advantages of 
simplicity and low-cost implementation. 
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[0013] According to another embodiment of the invention, the 
• light intensity distortion introduced by the system is also 
corrected by generating a light-intensity transfer function at 
each pixel and calculating coefficients for transfer arrays 
5 designed to normalize the system to produce a uniform response 
across each image and among images. This calibration procedure 
is carried out after geometric calibration and using measurement 
data already corrected for geometric distortion. The transform 
arrays are then applied to the measured intensities to generate 
10 uniform images. 

[0014] According to yet another embodiment of the invention 
relevant to interf erometric phase measurements, multiple phase 
measurements are taken and a random phase offset in the reference 
path length is introduced at each measurement. The multiple 
15 phase data so derived are then averaged to reduce phase-dependent 
systematic measurement errors. 

[0015] Other objects, features, and advantages of the invention 
will become apparent to those skilled in the art from a 
consideration of the following detailed description taken in 
20 conjunction with the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



[0016] FIG. 1 is a schematic view of a generic multi-channel 
imaging system. 

[0017] FIG. 2 is a schematic view of a multi-channel imaging system 
5 according to the prior art. 

[0018] FIG. 3. is a block diagram of the method of the present 
invention for calibrating and correcting for errors in a multi- 
channel imaging system. 

[0019] FIG. 4. is a diagram of the steps involved in determining 
10 image geometric distortion according to the invention. 

[0020] FIG. 5a is a schematic view of an exemplary apparatus for 

producing a calibration pattern. 

[0021] FIG. 5b. is a schematic view of an exemplary pattern and 

coordinates used to determine image distortion. 

15 . [0022] FIG. 6., is. a schematic view of a, surface used to visualize 
the distortion of each image. 

[0023] FIG. 7. is a diagram of the steps involved in correcting for 
image distortion. 
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[0024] FIG. 8. is a diagram of the steps involved in measuring 
intensity non-linearity in the sub-images. 

[0025] FIG. 9. is a diagram of the steps involved in correcting 
for image geometric distortion and intensity errors in a multi- 
channel imaging system. 

[0026] FIG. 10 is an. example of phase offsets and sampling periods 
used to reduce systematic measurement errors in interf erometric 
phase measurements. 
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DETAILED DESCRIPTION OF THE INVENTION 



[0027] The present invention provides apparatus and methodology 
for quantifying and correcting for errors in multi-channel imaging 
systems. There are three primary types of errors that are 
addressed in this invention: 1) geometric errors, such as image 
distortion and registration; 2) deterministic errors, such as 
channel or pixel dependent attenuation; and 3) stochastic or 
multivariate errors, such as thermal noise or sensor smear. For 
the purposes of this disclosure, the term "geometric distortion" is 
used to refer to those errors that produce distortion and 
misalignments in the images produced by the multi-channel system. 
The term "intensity distortion" is used to refer to nonunif ormity 
of intensity detection resulting from attenuation along the optical 
path, and detector nonlinearity or nonunif ormity of response. 

[0028] Turning to the drawings, a generic multi-channel imaging 
system is shown in FIG 1. The system consists of an image 
splitting element 1 that produces a plurality of image copies 20 
from an input image 10. The multi-channel imaging system can also 
be thought of as dividing an incoming wavefront 11 into a plurality 
of wavefront copies 14. The positions of the input image 10 and 
the image copies 20 are defined to be at conjugate image planes. 
The image copies 20 or wavefront copies 14 can be detected by a 
single or by multiple pixilated detector arrays. 
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[0029] FIG- 2 illustrates a prior art multi-channel imaging system 
used for measuring the phase difference between two substantially 
collinear, orthogonally polarized beams combined into an input 
wavefront 4. An input lens 12 focuses the orthogonally polarized 
5 beams onto a diffractive optical element 13 that produces a 
plurality of sub-wavef ronts 14, which are then further imaged by an 
output lens 15 onto a single detector array 17. A polarization 
phase plate 16 imparts different relative phase-shifts between the 
two orthogonally polarized beams and generates a plurality of 

10 phase-shifted interf erograms 6. Geometric image distortion, 
polarization rotation, and intensity non-uniformity can be 
introduced by the output lens 15, the polarization phase plate 16, 
and the detector array 17. As is well understood in the art, the 
common path elements, such as the input lens 12 and the diffractive 

15 element 13, can also introduce geometric distortion common to all 
images, as well as non-common distortion and intensity non- 
uniformity, if the orthogonally polarized beams deviate from being 
exactly collinear. 

[0030] The prior-art example shown in FIG.^ 2 is intended to 
20 exemplify how the image splitting element 1 of FIG. 1 is a 
combination of common path and non-common path elements that can 
introduce geometric distortion, polarization rotation, phase errors 
and intensity non-uniformity into the sub-images 20. The subject 
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of the present invention is a method for calibrating these types of 
errors and compensating for them during subsequent measurements. 

[0031] FIG. 3 illustrates the broad steps followed by the method of 
the invention for quantifying error sources in a multi-channel 
measurement. A first stage 100 consists of measuring field- 
dependent geometric distortion, followed by a second stage 110 
during which such field-dependent geometric distortion is applied 
to correct subsequent images, a third. stage 120 during which the 
light intensity transfer function is measured at each pixel, and 
finally a fourth stage 130 wherein the coefficients so generated 
are saved. An additional stage 140 is illustrated to indicate that 
the entire procedure is preferably repeated over a range of 
wavelengths . 

[0032] FIG. 4 illustrates the steps involved in stage 100 of the 
invention for quantifying and storing the image distortion present 
in the system. The first step 200 in the process is to project a 
pattern into the multi-channel imaging system. The pattern can be 
self-illuminated (e.g., light emitting diodes), back illuminated 
(e.g., light source and transparency), or front reflected (e.g., a 
light source and a reflective pattern) . A known pattern can be 
used to achieve absolute distortion correction, or a random, 
unknown pattern can be used to achieve relative alignment between 
the sub-images. For many real-world optical systems, the relative 
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alignment of images is the most significant correction needed to 
improve measurement accuracy because the absolute geometric 
distortion is acceptably small for comparison with the real world 
object. One example of this is found in interf erometric systems 
measuring strain on diffuse reflecting surfaces. The phase of an 
optical beam reflected from a diffuse surface varies rapidly with 
spatial position; therefore, the images must be accurately 
registered with respect to each other. However, the rapidly 
varying phase is subtracted out in subsequent measurements and only 
a slowly varying change in the surface height is measured. Because 
the measured quantity varies slowly with spatial position, the 
absolute geometrical distortion of the image may present an 
acceptably low error when comparing the results with other spatial 
data, such as a computer generated finite element model. 

[0033] In step 50, the calibration pattern is detected and 
digitized, for example by a CCD camera, a frame grabber, and a 
computer. The resulting images are then digitally parsed into N 
sub-images in step 60. The process of detection and parsing is 
described in the prior art. For example, in the prior-art example 
shown in FIG. 2, there are 4 sub-images. 

[0034] As shown in FIG. 4, the process followed for quantifying and 
storing the image distortion present in the system further includes 
steps 240, 260, and 280, respectively, for parsing the N sub-images 
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into j x k sub-regions where each sub-region has a center spatial 
coordinate given by x^,y k ; finding the absolute or relative offset 
of the pattern in each sub-region Ax N#j , k and Ay Nfj#k ; and performing 
a polynomial surface fit such that 

5 A*AU ( x j,k » yj,k ) = X C AU.».«*M and ( 1 ) 

Ay*.* ( x j,k » yj t k ) = 2 d NXnm x Jk y jk ( 2 j 

where c N , x#n , m and d N/X , m/n are coefficients of the smooth polynomial 
10 functions obtained in the fit. The subscript X denotes that the 
calibration may be performed at a plurality of wavelengths- The 
coefficients are then stored in step 300 and can be subsequently 
recalled and used for the correction of measurement data during the 
correction stage 110. The order of the polynomial functions and 
15 thus the number of coefficients can be adjusted in conventional 
manner as needed to obtain a good fit. 

[0035] FIG. 5a illustrates a calibration device 205 that can be 
used to produce a calibration image suitable for both absolute and 
relative image correction. The device consists of a substrate 220, 
20 which can be either transparent or coated with a reflective 
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coating, with a mask pattern 210 that could be a self-supported 
structure such as a thin metal or ceramic substrate patterned with 
repeated regions 212 of arbitrary shape. The mask pattern 210 can 
also be an integrated structure or coating such as an organic or 
5 inorganic thin-film deposited directly on the substrate 220 in such 
a way as to produce reflective or absorptive repeated regions 212 
. of arbitrary shape. For simplicity, only 9 repeated regions are 
. shown in FIG. 5a, but any number of repeated regions may be used as 
necessary to adequately characterize the multi-channel system. For 

10 typical commercial arrays of 1000 x 1000 pixels and a relative 
system distortion of about 0.5%, a 20 x 20 pattern can provide 
adequate sub-pixel resolution, although more points will typically 
lead to better correction. To produce the calibration pattern, the 
calibration device 205 is illuminated with a light source and 

15 either transmitted through the device or reflected from the 
surface. In either case, the illumination source should be 
arranged so that the light incident on the multi-channel imaging 
system is at substantially the same angle that will be used under 
actual measurement conditions. 

20 [0036] FIG. 5b illustrates how a sub-image 20 produced using the 
calibration device 205 of FIG. 5a is further analyzed according to 
the invention. The sub-image 20 comprises repeated features 26 
that may be either a bright spot surrounded by a dark background or 
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a dark spot surrounded by a bright background. The shape of the 
spot(s) may be important depending on the algorithm used, as 
described below. The sub-image is then . sub-divided into j by k 
sub-regions 22 such that each region contains one or many of the' 
repeated features 26. The location of the repeated features 26 
within each sub-region 22 is then determined using one of several 
methods described below. The location is defined by a set of x and 
y coordinates 28 and 29, respectively, that are referenced to a 
consistent origin point 25 within each sub-region 22. For example, 
the x and y coordinates 28 and 29 could be referenced to an origin 
point 25 located in the lower left-hand corner of each sub-region 
22, as illustrated in FIG. 5b. 

[0037] Several methods can be employed to determine the location of 
the repeated features 26. One is the use of a simple peak-finding 
routine to locate the brightest or darkest region in the sub-image. 
This requires that one repeated feature 26 be present in each sub- 
region 24. In a preferred embodiment of the invention, a centroid 
calculation is used, which is analogous to finding the center of 
mass within, the region . ^ This method .can ..produce a. more accurate 
result than the peak-finding method and it is permissible to have 
more than one repeated feature 26 within each sub-region 22. Each 
pixel's weight is the product of its relative x and y coordinates 
28 and 29 and the intensity value recorded at the pixel. The 
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centroid is then found by summing the products over all pixels and 
dividing the sum by the width of the region. Mathematically this 
can be expressed by, 

widthx 



kv - *J*re&» (3) and 



AY NX . k= H^l __ ', (4) 

N ' A ' M widthy 



where widthx and widthy are the width of each region along x and y, 
respectively, and I N/A (x,y) is the measured intensity value at each 
10 x/y pixel, in each sub-region N, and at each wavelength A. 

[0038] An additional step that can be used to further improve the 
accuracy of determining the location of the repeated features 26 
involves the use of a cross-correlation algorithm. In this case 
each sub-region is first cross-correlated with a master sub-region 
15 image according to the following formula: 



W( x '^ = 2I / ^ (ffl ' B) W('- m ^"" ) (5) 

ffl n 
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[0039] The master image can be taken from one of the sub images 
(relative correction) or can be synthesized from an ideal image 
(absolute correction) . The additional subscripts, i and j, for the 
intensity values denote the fact that the cross-correlation is 
5 performed separately for each sub-region of the image and not for 
the entire image at once* The resulting operation produces a 
unique correlation peak in each sub-region of each sub-image. 

[0040] A centroid or peak-finding algorithm can then be. applied to 
find the center for each sub-region of each sub-image. This method 
10 can produce excellent results even when there is significant noise 
present in the images and when there is significant distortion 
within the cell. 

[0041] Another approach suitable for the invention is to use a 
cross-correlation method with an unknown or random calibration test 

15 pattern. With this method it is only possible to determine the 
relative correction between images and not an. absolute calibration, 
but it has the advantage of working with a very simple, low cost 
target, for example a diffuse reflector or a ground glass 
transmission plate. One of the N sub-images is selected as the 

20 target or master image. Each sub-region from the master cell is 
cross-correlated with corresponding cells from the other sub-images 
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and a relative offset is determined. The other sub-images will be 
re-mapped to correspond to the chosen sub-image. This method is 
particularly useful when the overall distortion of the master image 
is low. FIG. 6 illustrates a 3-dimensional visualization of the 
5 error functions Ax N or Ay N for an arbitrary optical system. 

[0042] FIG. 7 shows the steps used in stage 110 to generate a 
correction map and apply it to subsequent data measurements. The 
stored coefficients are loaded in step 400 and transformation 
arrays are generated in a subsequent step 420 using the equations: 

10 

«,m 

15 [0043] The transform arrays can be pre-calculated to save 
processing time during the subsequent step 480. of correcting new 
data for distortion, which may be repeated multiple times, once for 
each new measurement. The distortion correction of step 480 
consists of interpolating a corrected data array from the raw input 

20 data array. A simple interpolation routine may be used, know as a 
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bi-linear transform, which consists of -finding the four nearest 
neighbor data points for each point in the original array and using 
a linear weighting to find the transformed value. Greater accuracy 
can be obtained using a higher order interpolation routine, such as 
a bi-cubic interpolation, which applies a piecewise polynomial, 
function to a 4 x 4 neighborhood of nearby points. 

[0044] The use of smoothly varying functions permits correction 
factors to be quickly calculated at any point in the field even 
though only several points were used to measure it. By selecting 
the polynomial order, the accuracy of the fit operation can be 
balanced against the spatial frequency response to mitigate the 
effects of random noise present in the measured data. For example, 
for many real-world cases a third-order fit produces an adequate 
balance between noise suppression and adequate spatial resolution. 

[0045] FIG. 8 illustrates the stage 120 of measuring the 
deterministic systematic errors. It comprises the step 500 of 
projecting a pattern, which could be either a uniform or "flat 
field" pattern or a spatially non-uniform pattern such as a 
sinusoidal intensity map. The incident pattern is digitized in 
step 50, parsed in step 60 and then subjected to the distortion 
correction algorithm of step 110 to resample the data. The image 
data are then processed by a threshold detection algorithm to 
identify bad pixels in step 520 and the intensity is used to fit 



coefficients of a single valued equation in step 540. An example 
of a simple linear fit is given by: 

I'„(x,y) = G N (x,y) I„(x,y) + 0 H (x,y), (8) 

where G N and 0 N are coefficients that quantify the relative gain and 
offset corrections necessary for each pixel. The coefficients are 
determined such that 

I'x(x f y) = I' 2 (x,y) = I' 3 (x,y) = I' 4 (x,y). (9) 

The input intensity is adjusted 560 and the loop is repeated 
until the coefficients converge to steady values. Once the 
values converge the intensity normalization coefficients and bad 
pixel maps are stored in step 600 for subsequent recall and use. 

[0046] FIG. 9 illustrates the approach of the invention to 
improve measurement accuracy after the calibration procedure is 
performed. The measurement is sensed by the detector array in a 
first step 40, digitized in step .50, parsed into N sub-images in 
step 60, resampled to correct for geometric distortion in step 
110, corrected for intensity non-linearity in step 700, corrected 
at bad pixels in step 800, and finally passed onto the actual 
measurement calculation in step 900. The actual process of image 
calibration, indicated by reference numeral 1000 in the figure, 
does not necessarily need to be performed in a sequential 
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fashion. That is, the steps cain be combined into, an efficient 
loop that only passes over each sub-image array a single time and 
performs the necessary calculations at each pixel. 

[0047] Bad pixels can also be corrected simply by setting a flag 
5 or code in the data to identify them as such and drop them from 
the calculation. However, for some systems it is possible to 
- reconstruct the bad or missing data from the corresponding pixels 
in the other sub-images. For example, in the case of four sub- 
images, this can be modeled by the relations: 

10 * I'i(x,y) = f 1 (I' 3 (x,y), I' 2 (x,y), I\(x,y)) 

I' 2 (x,y) ■= f 2 (I' 4 (x,y), I' 3 (x,y), I'^x^)) 

I' 3 (x,y) = f 3 .(I'i(x,y), I' 4 (x,y), I' 2 (x,y)) 

I' 4 (x,y) = f 4 (I' 2 (x,y), IS(x,y), I' 3 (x,y)) 

• where f n is a function that depends on the physical quantity that 
15 is measured. The advantage of reconstructing missing. data from 
good data in corresponding sub-images is. that the data can be 
processed with the downstream algorithm without passing of 
invalid data flags, which can slow processing or may not be 
handled by commercially available routines. In addition, data 
20 from all the sub-images can be displayed as contiguous. 

[0048] For example, for the case of simultaneous phase-shifting 
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interf erometry, where each sub-image is an interf erogram phase- 
shifted by 90 degrees, the functions f n are identical and given 
by: 

I' 1 (x / y) = I' 2 (x,y) + I f 4 .(x,y) - I' 3 (x,y) . 

5 I' 2 (x,y) = I'!(x f y) + I' 3 (x,y) - I' 4 (x,y) 

I' 3 (x,y) = I' 2 (x,y) + I' 4 (x,y) - I^x^y) 

I' 4 (x,y) = I'x^y) + r 3 (x,y) - I / 2 (x,y) 

[0049] Both geometric distortion and intensity distortion in a 
system may be s.trongly dependent on the illumination wavelength. 

10 These chromatic effects can be removed from the system by making 
calibration measurements at multiple wavelengths (as indicated by 
step 140 in FIG, 3) and subsequently recalling the proper 
calibration file as a function of operational wavelength for use 
for image calibration (see reference numeral 1000 in FIG. 9). 

15 The calibration coefficients can be interpolated to find a "best 
approximation" for a wavelength where a calibration procedure was 
not performed. 

[0050] Even after careful calibration, measurement errors will 
exist due to stochastic and multivariate error sources. For 
20 example, thermal noise in the electronics produces a random error 
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signal. According to yet another aspect of the invention, it is 
found that true random error sources can be reduced through 
averaging multiple measurements. Multivariate error sources are 
more problematic and, in general, cannot be reduced by averaging. 
One example of a multivariate error source that can cause 
problems when using a CCD camera is sensor smear, which results 
from light leakage into the vertical shift registers as a charge 
packet is clocked though a brightly illuminated region of the 
sensor. For interf erometric sensors, smear can cause a phase- 
dependent error signal that also . depends on the absolute 
intensity, fringe frequency and. orientation. Thus, it is not 
feasible to perform an a priori calibration for this type of 
error. However, because the error source manifests itself as a 
periodic function in a controllable variable, it is possible to 
induce a random nature to the error. 

[0051] In an interferometer, for example, the error can be 
randomized by shifting the phase offset between measurements. 
Averaging can then be used to reduce the error source, similar to 
stochastic noise. The measured phase of the interferometer can 
be expressed by 

(10) 
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where <p abs is the absolute phase difference between the reference 
and test beams, and e is the measurement error which depends 
explicitly on the absolute phase. The absolute phase difference 
can be split into an offset phase <p off that is constant over the 
surface and a spatially varying phase <p spa tiai according to 

so that 

tmeas (*> >0 = t spatial (*> >0 + Atf + ^spatial (*> >0 + ) . (12) 

The periodic error signal can be represented by an expansion 
series about the phase , 

^)=^a w cos(n^)+6 #f sin(w^) , (13) 

n 

where a n and b n are coefficients that may be generated on the 
basis of various factors such as smear, channel loss, etc. The 
subscript n corresponds to a harmonic of the fundamental 
frequency. 

[0052] For a static system, one without any phase drift or 
vibration, it is possible to introduce systematic phase delays 
and solve for the coefficients, a n and b n , using a least square 
fitting; however, measurements of practical interest for multi- 
channel systems are characterized by significant vibration and 



phase drift. It is also possible to fit the coefficients from a 
single inter ferogram provided sufficient phase periods (or 
fringes) are present (Schwider et. al.). However, the 
multivariate nature of the coefficients, that is, their 
5 dependence on multiple factors, such as intensity level and 
fringe orientation, makes it necessary to repeat the fitting 
process for every measurement. In addition, the process does not 
work for null-fringe conditions, which is of considerable 
practical importance. 

10 [0053] According to the invention, multiple measurements taken 
with a random phase offset between each measurement are averaged, 
as indicated by steps 900 and 1100 in Fig. 9. The final averaged 
phase at each surface coordinate x,y is given by 

15 where <p m is a random phase offset with uniform distribution 
between 0 and 2n. Both the phase offset and the periodic 
measurement error average to a constant value for all pixels. 
That is, 

j N 

TfY^^^^(x 9 y)+^J]^^^\£\^const, (15) 

20 where | e | is a constant value proportional to the magnitude of 
the periodic measurement error. Therefore, the resulting 
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averaged measurement is given by, 



(16) 



[0054] The constant value can easily be determined through spatial 
averaging of the data. The constant value can be subtracted to 
recover only the spatially varying phase, which is typically the. 
measurement of interest. . The periodic error is reduced by one over 
the square root of the sample size (l/\/N) , as is the stochastic 
noise. The method works for all fringe orientations and even in 
the presence of external periodic phase perturbations and phase 
drift. Therefore, it presents a robust solution for all potential 
situations. 

[0055] In order to introduce the phase offset, the reference or 
test optic can be mounted on a linear transducer to introduce a 
path-length delay equal to or greater than one wavelength. 
Equivalently, an electro-optic device can be inserted into the path 
of the reference or test beams to introduce the variable path 
length delay. Care must be taken to ensure that the phase offset 
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is uniform across all pixels. 



[0056] There are several ways to achieve a random phase offset 
distributed evenly over all phase values. In one embodiment the 
data sampling is set at a constant interval and a random phase 
5 offset is introduced temporally. A second, preferred, embodiment 
uses a random sampling interval while the phase offset is ramped 
with a constant period. This method is very easy to achieve in 
. practice since computers can easily generate an arbitrarily long 
sequence of random numbers with a uniform distribution to be used 
10 as a sampling delay. Similarly, a fixed frequency electronic ramp 
signal is also straightforward to produce. 

[0057] In a third embodiment, both the phase offset and sampling 
interval are periodic; however, the two intervals are adjusted to 
• be different from each other, and preferably non-integer multiples. 
15 The phase offset may be driven with sinusoidal excitation or, 
preferably, with a linear ramp. This method produces, quasi-random 
phase-offsets; however, it can produce an acceptable . approximation 
and adequate reduction in the error signal for some applications. 
The advantage of this embodiment is its simple implementation. 

20 [0058] An example of the phase offset and sampling periods for the 
preferred embodiment is shown in FIG. 10. The ramp function should 
be adjusted to introduce an optical path-length difference of (n 
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2n) , where n is an integer. The computer is used to generate a 
series of random delays having a uniform distribution over the ramp 
period. Other combinations of phase-shift waveforms and sampling 
frequencies can be employed in order to affect a uniform 
5 probability distribution between 0 and 2n phase offset. 

[0059] Those skilled in the art will understand that the 
preceding exemplary embodiments of the invention provide the 
foundation for numerous alternatives and modifications thereto. 
These other modifications are also within the scope of the 
10 invention, which is • accordingly not to be limited to the 
embodiments precisely shown and described above. 
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